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This article presents a modified output-prediction procedure and a new controller 
design based on the predictive control law. Also , a new predictive estimator is 
developed to complement the controller and to enhance system performance . The 
predictive controller is designed and applied to the tracking control of the Deep 
Space Network 70-m antennas. Simulation results show significant improvement in 
tracking performance over the linear-quadratic controller and estimator presently 
in use. 


I. Introduction 

The recent pointing requirements for the X-band 
(8.4 GHz) frequency, 70-m aperture antenna (further de- 
noted as the DSS 14 antenna) and the expectation of future 
K-band (32 GHz) capability dictate a need for high- 
performance controllers for the azimuth and elevation 
drives. This article presents a new design procedure for 
a tracking controller that significantly improves antenna 
tracking performance. It considers on-axis (or servo) 
tracking. In this case, the output is taken on the encoder 
(or tachometer) rather than the radio frequency (rf) point- 
ing position. The predictive controller uses future values 
of the stored input command to generate the control sig- 
nal. For this reason, the predictive control principle is 
considered useful in design of the Deep Space Network 
(DSN) antennas’ tracking controllers, since the antenna 
future tracking command is known when following stars 
or spacecraft. 

The tracking-control problem is a nontrivial exten- 
sion of the regulator problem, widely investigated in the 


control literature. Several approaches to the solution of 
the tracking problem have been presented in [1-11]. Pre- 
dictive controllers are described and analyzed in many pa- 
pers, among them [1,4,5,0,7,9,10,11]. In most of them, 
controlled auto- regressive and integrated moving-average 
(C A RIM A) models are developed and extensively used for 
output prediction and predictive control of linear systems 
[1,4,5,6,7,9,10]. State-space description serves as a stan- 
dard tool for system analysis and design. Besides tool stan- 
dardization, the state-space representation of a predictive 
control system provides a unique insight into system prop- 
erties, improves system design, and simplifies analysis. In- 
terpretations of CARIMA modeling in the state space are 
provided by [1,4,1 1]. 

This article presents a state-space description of the 
output-prediction procedure and a new controller design 
based on the predictive control law. A new input-reference 
scheme that uses the input-reference horizon is introduced. 
Thus, the increment of the control signal is determined 
with respect to the input horizon rather than to the last 



value of the input. Also, the introduced weighting matrix 
includes a forgetting factor. Both features significantly 
improve system performance. 

The usefulness of the predictive control law depends on 
the availability of the plant-state variables for measure- 
ment. Typically, not all state variables can be measured, 
although for an observable system they can be estimated. 
In this article, a predictive estimator is developed as a 
complement to the predictive controller to speed the esti- 
mation process and enhance system performance. 

Using the introduced predictive control and estimation 
laws, the state-space predictive controller is designed and 
applied to the tracking control of the DSN 70-m antennas. 
Simulation results show significant improvement in track- 
ing performance over the linear-quadratic (LQ) controller 
and estimator presently in use. The robustness to param- 
eter variation and the disturbance-suppression properties 
are found to be fairly good for the considered predictive 
control system. 


II. Output Prediction for a Linear System 

A plant with nu inputs and ny outputs is considered. 
Its linear state-space model consists of n states 


first NR time instants and is constant for the remaining 
NU — NR instants 


U r (k) 


ii t -_ i(Jk) for Jfc = \ } ...,NR 
Ui-i(NR) for jfc = NR + U . . . , NU 


(3) 


where the integer N R < NU is the lengt h of the reference 
horizon. Thus, 


Ur(i) = EU(i-l) (4) 
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NR times 
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>NR times 


and I nu is the identity matrix of dimension nu. 


x(i -f 1) = Ax(i) -f Bu(i), y(i) = Cx(i) + Du(i) (1) 

where x E R n , u E R nu , and y E R ny . The task is to pre- 
dict output y for NY steps ahead, given projected input 
u for NU steps ahead. The integer NY is the length of 
the output horizon, while NU is the length of the input 
horizon. For casual systems, the length of the input hori- 
zon does not exceed the length of the output horizon, i.e., 
NU<NY. 

Before determining the predicted output, the input and 
output sequences (or horizons) are introduced. Three 
types of input sequences are defined. First is the input 
horizon U(i), consisting of the input from an instant i up 
to NU — 1 steps ahead 

U T (i)= [uf(0),nf(l),...,uJ(NU -1)] (2) 

where U((k) is the predicted input at instant i with k steps 
ahead. The input horizon U(i — 1) is a horizon predicted 
at the previous time instant. Note that it is not a delayed 
prediction at instant i, i.e., Ui_i(k) / Ui(k — 1). 

The second sequence, reference-input horizon U r (i) } is 
identical to the previous input horizon U(i — 1) for the 


The last input sequence, input-increment horizon 
A£/(i), is defined with respect to reference horizon U r (i) 
as follows: 

AU(i) = U(i) - U r (i) = U(i) - EU(i - 1) (6) 

Sequences U(i), U(i- 1), U r (i), and A U(i) are shown in 
Fig. 1. 

Two output sequences are introduced: output horizon 

y 

and predicted-output horizon Y 

Y T (i) = lyI(l),yJ(2),.-.,yI(NY)} (8) 

The latter is an output of the system with reference hori- 
zon U r as an input. Components */,(&) and y t (i b) are the 
output and predicted output, respectively, at instant i with 
k steps ahead. Note that although the system output at in- 
stant i with k steps ahead is equal to the output at instant 
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i -f l with k - / steps ahead, y;(fc) = Vi+i(k - 0 = J/(* + 0* 
/ < Ar , the same is not true for the predicted output. The 
prediction at instant i with k steps ahead is not the same 
as the prediction at instant i + / with k — l steps ahead. 

The output horizon is obtained from the plant model, 
Eq. (1), for k = 1, . . . , NY: 

y(i + k) = CA k x(i) + CA k ~ l Bu(i) + ... 

+ CBu(i + k - 1) (9) 

Predicted output Vi(k) is defined as a system response 
to the reference-horizon input U r (i )• Thus, for k = 

y € (Jb) = + + (10) 

and Eq. (9) is now 

»(i + fc) = y(fc) + CA^'BAuiiO) + ... 

+ (CA k - NR ~ 1 B + .. . + CB)Aui(NR- 1) 

( 11 ) 

Denoting the Markovian matrix 

[9i 0 0 1 


^ I 

92 

9 1 

0 



( 12 ) 

9nu 

9 NU — 1 

9l 

9ny 

9NY - 1 

■ ■ • 9NY-NU + 1 - 



where y, = CA X 1 B is the ith Markov parameter, one 
obtains the output horizon from Eq. (11): 

Y(i) = Y(i) + GAU(t) (13) 

Predicted-output horizon Y, necessary to determine Y in 
Eq. (13), is determined from Eq. (10): 

Y(i) = Hx(i) + GU r {i) = Hx{i) + GEU(i - 1) 

= Hx(i) + FU(i- 1) (14) 


where 

F = GE y II = O NY -iA (15) 

E is as given by Eq. (5), and 

Osr-i = [C T (CAf...(CA NY -r] T (16) 

The input-increment horizon depends on the length of 
the reference horizon. In particular, for NR = 1, one 
obtains Aui(/:) = Ui(k) — U{ _ i(l), k = 1, . . . , NU. This 
is the case of the generalized predictive control of [5], where 
the control increments are defined with respect to the last 
input command. For NR = NU y one obtains E = I 
and AC/(i) = C/(i) — £/(i — 1). In this case, the control 
increment is defined with respect to the previous control 
over the whole length of the input horizon. If the input 
increment is determined with respect to the zero-reference 
input, the input-increment horizon is equal to the input 
horizon, A U(i) = t/(z), and for this case, NR = 0. Signal 
sequences U(i) } U(i — 1), U r (i), and A U(i) are shown in 
Fig. 1 for NR =0,1, and NU and for a generic NR. Note 
also that for NR — NU one obtains F = G, for NR = 1 
one obtains F T = [gj , gj + gf , ... +^y]> and for 

NR = 0 one obtains F = 0. In the latter case, the output 
is predicted from the system state only, while otherwise it 
is predicted from the state and the system input as well. 

III. Predictive Control 

The basic task for a predictive controller is to assure 
that for the bounded input the future output Y will closely 
follow the input command Y 0 within the output horizon 
NY: 

Yj (0 = b£(l). y T oi{ 2). • • ■ , vh(NY)] (17) 

where y 0 i(k) is the command signal at instant i with k 
steps ahead. Thus, the task is to minimize the plant track- 
ing error while the input remains bounded: 

e(i) = n«-y(0 (is) 

e T (i) = [£f(l),ef(2),...,ef(A r y)] (19) 

where Si(k) is the error at instant i with k steps ahead. 
The tracking error within the output horizon, as well as 


32 



the restrictions on the input within the input horizon, are 
included in performance index J 

J =tr ( e T (i)Qe(i ) + AU T (i) RAU («)) (20) 

where tr(.) denotes the trace of a square matrix; Q and 
R are symmetric, positive, definite matrices; Q is the 
tracking-error weighting matrix; and R is the control-effort 
weighting matrix. 

The necessary condition for the optimum, dJ fdAU — 
0, applied to Eq. (20) and using Eq. (13), yields 

A U(i) = K(Y 0 (i) - 7 ( 0 ) = ia(i) (21) 

where 

I< = (G t QG+ R)~ 1 G t Q (22) 

and ?(*) = Y 0 (i) — Y(i) is the predicted output error. The 
resulting control increment AU(i) covers the whole input- 
horizon NU] for control purposes, however, only the first 
component (the current control increment) is used. Let k 
denote the first nu rows of I\ : 

k = eK , e = [7 nu 0 0 . . . 0) (23) 

Then the control increment at instant i is 

Au(t) = k(Y 0 (i) - 7 ( 0 ) = ke(i) (24) 

and the control input is obtained from Eq. (6) 

u(0 = u(i - 1) + Au(i) (25) 

Combining Eqs. (24), (25), (14), (15), and (16), one ob- 
tains the control command at moment i: 

u(i) = u(i - i) + * (Y 0 (i) - 7(0) 

= u(i - 1) -j- kY 0 (i) - kHx(i) - kFU(i - 1) 
and with u(i — 1) = eU(i — 1), the above equation yields 

ti(0 = *n(0 - *^*(0 + (C - kF)U(i - 1) (26) 


Thus, the command u(0 depends on the previous input 
horizon U(i — 1), on the actual state x(i), and on the 
control command Y 0 (i) up to NY steps ahead. 

The closed-loop system equations are obtained by com- 
bining the plant equation, Eq. (1), with the controller 
equation, Eq. (26), and by introducing the new state vari- 
able U 0 (i) such that 

U 0 (i+1) = U(i) (27) 

In this way, one obtains 

x(i + 1) = (A - BkH)x(i) + B(e - kF)U 0 (i) + BkY 0 (i) 

U 0 (i + 1) = -KHx(i) + (I N - I<F)U 0 (i) -f I<Y 0 (i) 

y(i) = Cx(i) (28) 

and N = NU x nu. With the new state variable z T = 
[a: T ,Uj], the closed-loop equations are 

z(t + 1) = A c z(i) + B c Y 0 (i), y(i) = C c z(i) (29) 

where 

. _ \A- BkH B(e-kF)' 

Ac ~[ -KH I n -ICFF 


_\l n -BKO NY -i)A B(e-kF) 1 
-KH I N - KF J 

Bj = [ Bk I<) T , C c = [C 0] (31) 

One can see that control command ti(i’) is now 

u(i) = kY 0 (i) + [~kH e - kF]z(i) (32) 

fully recovered from the current state of the system and 
from the input command. 

The block diagram for the closed-loop system, 
Eqs. (29)— (3 1) , is presented in Fig. 2. The system consists 
of the plant, the predictor (PRD), the controller (CO), 
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and the command horizon generator (CHG). The predic- 
tor structure is shown in Fig. 3(a) and the command hori- 
zon generator in Fig. 3(b). The state-space representation 
of the command horizon generator is 


variables) and input disturbances i> u (i) (Fig. 2). The dis- 
turbances are included in the closed-loop system model, 
with the triple ( A Cy B yy C c ) for the output noise and the 
triple (A c , B u , C c ) for the input disturbances, where 
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where I is an identity matrix of order ny, and A r h and B r h 
have NY - 1 rows, while C r h and D r h have NY columns. 


Their impact on system performance is studied in Sec- 
tion V.D. For high-frequency disturbances, the distur- 
bance-rejection properties of the system significantly im- 
prove when the lowpass filter (LPF) is applied as in 
Fig. 4(a). The plant states related to the command sig- 
nal are obtained from the plant model (PM), and they are 
extracted from the measured states. The resulting signal 
passes through a lowpass filter and is added to the states 
previously extracted. The filter is shown in Fig. 4(b). 


The weighting matrices R and Q are the tuning param- 
eters of the optimal design (for example, see [2,12]). That 
means they are to be adjusted until satisfactory results 
are obtained. Although they are not “active” in the op- 
timal design solution, their choice significantly influences 
the performance and stability of the system. A general 
procedure for a reasonable choice of the weighting matri- 
ces is not yet known. In this article, a simplified procedure 
is developed. The weighting matrices obtained from this 
procedure significantly improve system performance, i.e., 
tracking error. 

A diagonal matrix R = pi has been chosen as a control 
weighting matrix, where p > 0 is a scalar, and tracking- 
error weighting matrix Q has the following structure: 

Q = diag(q,otq,a 2 q, . . . ,a NY ~ 1 q) (34) 

The diagonal component qa k ~ l q is the weight of the error 
of Ci(k) y the Arth component of £r(i). The last weight is time 
dependent; the weight of the output error at the (i + k ) th 
time instant is a k ~ l . The scalar a is called a forgetting 
factor. The most recent output is given a unit weight, and 
the future output penalized (in fact, awarded, as will be 
shown later) exponentially. With this arrangement, the 
choice of R and Q reduces to the choice of parameters p, 
a, and q, as is illustrated in Section V.B. 

There are two sources of system disturbances; mea- 
surement noise u y (i) (or v x (i) when measuring all state 


IV. Predictive Estimation 

Implementation of the predictive controller depends on 
the availability of the plant states for measurement. Often, 
these parameters are not available. An LQ estimator (for 
example, see [2,13]) that estimates plant state from its 
output can be considered as a solution to the problem. 
Its action, however, is too slow for the predictive control 
system, and the predictive scheme is included in the design 
of the estimator. Thus, a new estimator with dynamic 
characteristics comparable to the predictive controller is 
developed. 

The estimate x (i) of the plant state ar(i) is determined 
from the input and output horizons as follows. From 
Eq. (1), one obtains 

k 

CA k x(i) = y(i + k)-J2 CA’-'Bu(i +k- j), 

j = 1 

£ = 0,1,2,..., NY 

or 

0 N yx(i) = Y e (i)-G € U(i) (35) 

where 
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and H , G, and Y are given in Eqs. (15), (12), and (7), 
respectively. Variable Y e is an augmented output hori- 
zon composed of the current output y(t’) and the output 
horizon Y(i). From Eq. (35), the estimate £(i) of x(i) is 
determined such that for a symmetric, positive, weighting 
matrix Q e the estimation index 

J e - \\H e x(i) - H e x(i)\\Q e (37) 

is minimal, obtaining 

x(i) = 0+ y (Y e (i)-G e U(i)) (38) 

where O^y = 0 0 T NY Q e O NY )~ l 0 T NY Q e - Note that the 
state estimate is determined from input and output hori- 
zons, while input and output signals, rather than horizons, 
are available for estimation. The input horizon is available, 
however, right after the controller output (Fig. 2). Output 
horizon Y(i) is not available directly; nevertheless, it can 
be obtained from the plant model as follows: 

X(i 4 1) = AX{i) 4 BU(i), Y(i) = CX(i) 4 DU(i) 

(39) 

The estimator is shown in Fig. 5(a). Thus, the plant 
state is estimated from its output and the input horizon. 
This scheme is similar to the LQ estimation scheme, since 
it uses the available input and output signals and the plant 
model to generate the estimate. The block diagram of 
the predictive control system with the predictive estimator 
(EST) is shown in Fig. 6. 

Unlike the LQ estimator, the predictive estimator does 
not have filtering properties, since its output x(i) is pro- 
portional to a noisy signal y(i). This drawback can be 
removed as follows. Given the plant model output y n (i) 
the output error e y (i) = y(i) — y n (0 is filtered by a proper 
filter, obtaining the filtered error £ y j(i ). In most cases, the 
output error is a high-frequency noise; hence, a lowpass fil- 
ter is applied. The filtered output is obtained by adding a 
filtered error to the nominal output yj(i) = !/n(i) + £ y/( J T 
In this way, most of the noise power is removed from the 
output signal, while the basic properties of the signal re- 
main untouched. The estimator with a filter is shown in 
Fig. 5(b). The filter action will be illustrated in the next 
section. 


V. Predictive Control and Estimation for the 
DSS 14 Antenna 

Performances of the predictive controller and estimator 
are checked through tracking simulations of the DSS 14 
National Aeronautics and Space Administration (NASA)/ 
Jet Propulsion Laboratory (JPL) 70-meter DSN antennas. 
The existing control scheme for the DSN 70-m antennas 
[14] is based on an LQ regulator design with the integral 
action as presented in [15-20]. The LQ control system is 
shown in Fig. 7, in which the plant output is augmented 
by the addition of the output integrals in order to ensure 
the zero mean value of the constant-rate tracking error. 
The LQ controller is designed for this augmented plant 
with a constant tracking command. This assumption can 
be a significant source of tracking error. A controller de- 
signed for the constant tracking command can result in 
insufficient antenna performance, especially for relatively 
fast commands or varying rate commands. In this section, 
the performance of the predictive controller is compared 
with that of the LQ controller in the tracking environment. 

A. Plant Model 

The state-space model of the DSS 14 antenna [14] is a 
four-state model with position rate u as an input and po- 
sition rate y as an output. Its discrete-time representation 
(Ad> Bd,C<t), with the sampling period At = 0.05 sec, is 
obtained from the continuous-time representation in [2]: 

*0.0468 0 0 0 “ 

0 0.5443 0.3474 0 

Ad ~ 0 -0.3474 0.5443 0 

0 0 0 0.8872. 

-0.0113- 

B d = 0 0399 ’ ° d = [0 ' 7239 9 2260 ° 1H21] 

.0.0538. 

The system (Ad, B dl Cd) is augmented. As a result, its 
output consists of the position rate, the angular position, 
and the integral of the position. The augmented system 
is shown in Fig. 8. Denoting Xd the state of the system 
(Ad, Bd, C d ), and x po and x ipo the position and the integral 
of the position, respectively, one obtains from Fig. 8 

Xi p o(i T 1) = Xi po (i ) T At x po (i) 

x p0 (i 4- 1) = Xp 0 (i) 4 AtCdXd(i) 

Xd(i 4 1) = AdXd(i) 4 B d n(i) 
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With the plant-state variable x T = [xi p0} x po , arj], the 
triple (A,B,C) is the resulting plant-state space repre- 
sentation 
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used in simulations presented below. 

B. Weighting Matrices and Input and Output Horizons 

For simulation purposes, a piecewise-linear profile of the 
position command is chosen, with linear increase followed 
by linear decrease and the final constant value (Fig. 9). 
The command rate is 4 mdeg/sec, which is a typical si- 
dereal tracking rate. The shape of the command is more 
dramatic than the real tracking command, but it has been 
chosen to emphasize the tracking possibilities of the pre- 
dictive controller. A more realistic tracking command will 
be used later in this article. 

In order to perform a series of simulations, weighting 
matrices R and Q are chosen such that the output er- 
ror is small while the control effort is maintained within 
reasonable limits. For diagonal-control weighting matrix 
R = pi , parameter p = 0.01 is chosen. Tracking-error 
weighting matrix Q is as in Eq. (34). Component q is in 
the form q — diag(<?j, q p , q r ). It represents the weight of 
the integral, position, and rate components of the output. 
The following choices of weight are recommended from a 
series of simulations tracking the command as in Fig. 9: for 
the integral-of-the-position signal, g, = 10; for the position 
signal, q p = 1; and for the rate signal, q r = 0.1. Coeffi- 
cient a*” 1 in the weighting matrix is the weight of the 
k th error component in the output horizon. Simulations 
have been performed in order to determine the value of 
parameter a. The plot of the Euclidean norm of tracking 
error ||y — y 0 || 2 versus a is shown in Fig. 10 both for differ- 
ent lengths of output horizon and for lengths of input and 
reference horizon equal to lengths of output horizon. The 
plot shows the minimal tracking error obtained for a — 6.2 
and NY = NU = NR = 6. For NR = NU = NY, Fig. 11 
plots the values of a for which the tracking error is min- 
imal. The figure shows that, for horizons that are long 
enough, the forgetting factor is close to 1. Thus, in this 
case, the time weighting does not significantly improve the 
tracking error. However, for short horizons, the proper 
choice of forgetting factor is a critical factor that mini- 
mizes the error dramatically. In prediction, the forgetting 
factor is greater than one (thus, “reminding factor” could 
be an adequate name for it). This is in contrast to the 
forgetting factor value in estimation procedures such as in 


[12], where the factor is smaller than one. This difference 
occurs due to opposite time directions; in estimation, the 
past values of signal are weighted, while in prediction the 
future values of signal are processed. 

From simulations, the impact of the length of input- 
reference horizon NR and output horizon NY on the 
tracking error is determined. The results are plotted in 
Fig. 12. One can see that for NR > n/2 and NY > 2 n 
(where n = 6 is the number of plant-state variables), the 
performance error is close to the minimal one. 

C. Antenna Performance 

The performance of the DSS 14 antenna with the track- 
ing command as shown in Fig. 9 has been evaluated for 
the parameters recommended above. The following pa- 
rameters for the predictive controller were chosen: NR = 
NU = NY = n = 6 and weighting matrices with p — 0.01, 
q = diag(10, 1,0.1), and a = 6.2. The reference signal 
and the position of the antenna with the predictive con- 
troller for ct — 6.2, for a = 1, and for the antenna with 
the LQ controller are shown in Fig. 9. The prediction er- 
rors and control input for the above three cases (a = 6.2, 
a = 1, and LQ controller) are shown in Figs. 13 and 14. 
The figures show better performance by the predictive con- 
trollers than by the proposed LQ controller with compara- 
ble control effort. Also, predictive-controller performance 
with time-weighted output error (a > 1) is better than 
predictive-controller performance without time weighting 

= !)• 

The minimum of tracking error for output horizon 
NY = 4 is obtained for a = 6.4. These two parameters are 
used in further simulations, since it is reasonable to have 
the length of the output horizon as small as possible; the 
dimension of the controller as well as the complexity of the 
system depend on NY . The step-response and frequency- 
response plots of the closed-loop system with a predictive 
controller and an LQ controller are compared in Figs. 15 
and 16. Figure 15 shows that the settling time and over- 
shoot for the system with the predictive controller, with 
NY = NU — NR = 4 and a — 6.4, are significantly re- 
duced from the system with the LQ controller. Similarly, 
from Fig. 16, one can see the tracking performance is im- 
proved; the magnitude of the closed-loop transfer function 
is equal to 1 over a wider bandwidth. Also, roll-off rate is 
improved for the system with a predictive controller (with 
NY — NU = NR = 4 and a = 6.4, as well as with 
NY — NU — NR — 6 and a = 5) when compared to the 
system with an LQ controller. 

The piecewise constant-rate command, as well as the 
unit step command, are dramatic scenarios for the DSS 14 
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antenna and have been introduced in order to present dy- 
namic possibilities of the predictive control. In order to 
meet the typical working requirements for the antenna, 
the raised-cosine command is introduced, as in Fig. 17 
(solid line). This kind of command is close to the real ele- 
vation or azimuth trajectory of the antenna (conscan-like 
tracking). The plot of the output of the predictive con- 
trol system overlaps the plot of the command, while the 
output of the LQ control system is plotted by a dashed 
line (Fig. 17). The tracking error, the difference between 
the output and the command, is plotted in Fig. 18(a) for 
the LQ control system and in Fig. 18(b) for the predictive 
control system. For the LQ control system, the error is on 
the order of 10~ 4 , while the error for the predictive control 
system is on the order of 10“ 7 . In both cases, however, the 
control effort is almost the same (Fig. 19). 

D. Robustness and Disturbance Suppression 

The robustness of the closed-loop system to the plant- 
parameter variations is checked as follows. The plant poles 
are randomly perturbed within 20 percent margin, and 
the error in the step-command tracking is simulated for 
500 random samples. The results of the simulations are 
presented in Fig. 20. On the average, tracking error has 
changed about 5 percent in comparison with the nominal 
plant error, and the maximal tracking error is 66 percent 
larger than the nominal plant error. The step and fre- 
quency responses of the closed-loop system for the nomi- 
nal plant and for the plant model deviated 20 percent from 
the nominal, as is shown in Figs. 21 and 22. Both plots 
show good performance and robustness of the system. 

Two sources of disturbances of the antenna are stud- 
ied: the input disturbances and output disturbances (mea- 
surement noise). The input-disturbance transfer functions 
(from v u to y) and output-disturbance transfer functions 
(from v x to y) are shown in Fig. 23, the latter one for both 
position- and rate-measurement noise. One can see from 
Fig. 23 that the input disturbances are significantly sup- 
pressed, while the position-measurement noise is amplified 
over certain frequency ranges. 

The nature of the antenna disturbances is not sat- 
isfactorily known, and here their general properties are 
outlined. Input disturbances, such as wind or thermal 
forces, are low-frequency signals. Measurement noise, on 
the other hand, is a high-frequency signal (high in compar- 
ison to the antenna fundamental frequency, which is less 
than 1 Hz). Therefore, for testing purposes, white-noise 
input disturbances and high-frequency measurement noise 
with frequency components over 3 IIz are applied. The sys- 
tem response due to different signal-to-noise (S/N) ratios 


is simulated. The results are compared in Fig. 24, where 
good disturbance-suppression properties of the system are 
observed. The impact of the input noise is much smaller 
than that of measurement noise. This feature can be ex- 
plained with the lowpass- filtering property of the plant; 
the noise is filtered before entering the predictor. The 
tracking error and plant input for input noise with S/N 
ratio = 10 and for measurement noise with S/N ratio = 
100 are shown in Figs. 25 and 26. The effect of the mea- 
surement noise is reduced by applying a filter, as in Fig. 4. 
The transfer-function plots from the output disturbances 
to the system output for the system with the filter are 
shown in Fig. 27. The tracking error due to measurement 
noise is reduced significantly, even for white noise (Fig. 28). 

E. Predictive Estimator 

Predictive-estimator performance is compared to the 
performance of the LQ estimator. The plant model 
(Ad, Bd j Cd) has been used for simulations, with unit- 
step input and zero-initial conditions. For estimation 
purposes, the initial conditions have been changed to 
[0.1 0.1 0.1 0.1] T . The estimation results are shown in 
Fig. 29. The LQ estimator needs approximately 2 sec- 
onds to reach an acceptable estimation error, while the 
predictive estimator determines the states in virtually no 
time. In the case of noisy output, with S/N ratio = 100, 
one obtains estimation errors for the LQ and predictive 
estimators (with and without filter) as in Fig. 30. The 
unfiltered predictive estimate in Fig. 30(b) would make 
the estimator useless for prediction purposes. However, 
these errors are reduced by a filter, as in Fig. 30(c). The 
maximum error of the predictive estimator with a filter is 
much smaller than the residual error of the LQ estimator 
even after 4 seconds in action. Finally, simulations indi- 
cate that the identity-weighting matrix ( Q e = I) is the 
optimal choice for estimation purposes. 

VI. Conclusions 

In this article, a modified state-space predictive con- 
troller is introduced, and a predictive estimator is pre- 
sented to complement the design of a predictive-control 
law. This approach has been used for the design of the 
tracking controllers of the NASA/JPL 70-m antennas. 
Several tracking scenarios have been introduced (step in- 
put, constant-rate rise and fall, raised-cosine trajectory) 
to test the tracking behavior of the predictive controller. 
Significant improvement of performance for presented sce- 
narios has been observed. It has been shown that the time 
for the predictive estimator to reach an acceptable level 
of estimation error is much smaller than that for the LQ 
estimator. Also, a wider bandwidth and improved roll-off 
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rate is obtained for the predictive closed-loop system in 
comparison with the LQ regulator system. The predictive 
control system is robust to the plant-parameter variations. 
Shifts of plant poles of 20 percent of their nominal values 
keep the tracking performance good; the tracking error is 
on the same order as for a nominal plant. Disturbance- 


suppression properties of a predictive control system also 
have been simulated and found to be good for input distur- 
bances and measurement noise if the measurement-noise 
spectrum is higher than the plant-fundamental frequency. 
The system disturbance-suppression properties can be en- 
hanced if the disturbance filter is included in the system. 
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Fig. 8. The augmented model of the DSS 14 antenna. 



Fig. 9. The command signal and antenna output. 



a 

Fig. 10. Tracking error versus forgetting factor a. for different 
lengths of output horizon NY. 



Fig. 11. Forgetting factor a, for which the minimal tracking 
error Is achieved versus length of output horizon NY. 
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ig. 12. Tracking error versus: (a) length of reference horizon 
NR, and (b) length of output horizon NY. 



Fig. 13. Tracking error. 
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Fig. 16. Closed-loop system frequency response: (a) magnitude, 
and (b) phase. 
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Fig. 19. Antenna input. 
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Fig. 20. Tracking-error ratio of the 20-percent randomly deviated 
model to the nominal model. 
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ig. 21. Closed-loop step response with the 20-percent-deviated 
model. 
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Fig. 22. Closed-loop frequency response: (a) magnitude, and 
(b) phase. 
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Fig. 23. Transfer functions: (a) magnitude, and (b) phase. 




Fig. 25. White noise input disturbances with nolse-to-signal 
ratio = 0.1: (a) tracking error, and (b) antenna input. 


Fig. 24. Tracking errors due to input disturbances. 
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Fig. 28. Tracking error for white-noise c 
system with a noise 



Fig. 29. Estimation error: (a) LQ estima 
estimator. 









